Galactic Collapse of Scalar Field Dark Matter 
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We present a scenario for galaxy formation based on the hypothesis of scalar field dark matter. 
We interpret galaxy formation through the collapse of a scalar field fluctuation. We find that a 
cosh potential for the self-interaction of the scalar field provides a reasonable scenario for galactic 
formation, which is in agreement with cosmological observations and phenomenological studies in 
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In the last years, the quest concerning the nature of 
the dark matter in the Universe has received much at- 
tention and has become of great importance for under- 
standing the structure formation in the Universe. Some 
candidates for dark matter have been discarded and some 
others have recently appeared. The standard candidates 
of the Cold Dark Matter (CDM) model are axions and 
WIMP'S (Weakly Interacting Massive Particles), which 
are themselves not free of problems. Axions are mas- 
sive scalar particles with no self interaction. In order 
for axions to be an essential component of the dark 
matter content of the Universe, their mass should be 
m ~ 10 _5 eU. With this axion mass, the scalar field col- 
lapses forming compact objects with masses of order of 
M crit ~ 0.6mp;/m ~ 1CU 6 M fl],||], which corresponds 
to objects with the mass of a planet. Since the dark mat- 
ter mass in galaxies is ten times higher than the luminous 
matter, we would need tenths of millions of such objects 
around the solar system, which is clearly not the case. On 
the other hand, there are many viable particles with nice 
features in super-symmetric theories, such as WIMP'S. 
However, since these candidates behave just like stan- 
dard CDM, they can not explain the observed scarcity 
of dwarf galaxies and the smoothness of the galactic-core 
matter densities, since high resolution numerical simu- 
lations with standard CDM predict an excess of dwarf 
galaxies and density profiles with cusps Q. This is the 
reason why we need to look for alternative candidates 
that can explain both the structure formation at cosmo- 
logical level, the observed amount of dwarf galaxies, and 
the dark matter density profile in the core of galaxies. 

In a recent series of papers, we have proposed that the 
dark matter in the Universe is of a scalar field nature with 



a strong self-interaction The scalar field has been 

proposed as a viable candidate, since it mimics standard 
CDM above galactic scales very well, reproducing most of 
the features of the standard Lambda Cold Dark Matter 
(ACDM) model Pf^,P,|lCfl . However, at galactic scales, 



the scalar field model presents some advantages over the 
standard ACDM model. For example, it can explain the 
observed scarcity of dwarf galaxies since it produces a 
sharp cut-off in the Mass Power Spectrum. Also, its self- 
interaction can, in principle, explain the smoothness of 
the energy density profile in the core of galaxies J7|,pd|. 
Nevertheless, the main problem when a new dark matter 
candidate is proposed is the study of the final object that 
would be formed as a result of a gravitational collapse. 

The formation of galaxies through gravitational col- 
lapse of dark matter is not an easy problem to under- 
stand. A good model for galaxy formation has to take 
into account all the observed features of real galaxies. 
For example, it seems that many disc galaxies contain a 
black hole in their center, but others do not jl2). Typical 
galaxies are spiral, elliptical or dwarf galaxies (irregular 
galaxies may be galaxies still evolving). In most spiral 
and elliptical galaxies the luminous matter extends to 
~ 10 — 30kpc, and the total content of matter (includ- 
ing dark matter) is of the order of 10 10 — 10 12 M Q , with 
about 10 times more dark matter than luminous one. 
The central density profile of the dark matter in galaxies 
should not be cusp (l]| . Even though the luminous mat- 
ter represents only a small fraction of the total amount 
of matter in galaxies, it plays an important role in galaxy 
formation and stability [Q. On the other hand, it is still 
not well established if the mass of the central black hole 
and the mass of the halo are correlated |i4|] , etc. 

There are some ideas in this respect when dealing with 
a scalar field. It is known that the final stage of a col- 
lapsed scalar field could be a massive object made of 
scalar field particles in quantum coherent states, like bo- 
son stars (for a complex scalar field) or oscillatons (for 
a real scalar field) 0,^,|l5|. It is thus important to in- 
vestigate whether the scalar field would collapse to form 
structures of the size of galaxies and provide the cor- 
rect properties of any galactic dark matter candidate, 
like growing rotation curves and appropriate dark mat- 
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ter distribution functions. Also, we need to know which 
are the conditions that must be imposed on the scalar 
field particles. 

Our main aim in this paper is to present a plausi- 
ble scenario for galaxy formation under the scalar field 
dark matter (SFDM) hypothesis. Through a gravita- 
tional cooling process , a cosmological fluctuation of 
the scalar field collapses to form a compact oscillaton by 
ejecting part of the field. The key idea consists precisely 
in assuming that such final object could distribute as 
galactic dark matter does. The final configuration then 
should consist of a central object (a core), i.e. an oscil- 
laton, surrounded by a diffuse cloud of scalar field, both 
formed at the same time due to the same collapse process. 

At the cosmological scale, it is found that the mass of 
the boson is not the only parameter that determines the 
power spectrum. The self-interaction of the scalar field 
is also important. Following an analogous procedure to 
the one used in particle physics, we may write a phe- 
nomenological Lagrangian with all the terms we need in 
order to reproduce the observed Universe. In particular, 
if one uses a minimally coupled real scalar field $ with a 
self-interaction potential of the form (kq = 8ttG, we use 
natural units such that Ti — c — 1) 



l/($) = Vb[cosh(Ay^)-l] , 



(1) 



then one can show |g|7| that the SFDM model reproduces 
well all the successes of the standard ACDM above galac- 
tic scales. The free parameters of the scalar potential, Vq 
and the scalar field mass m$ = Aa/Vo^o, can be fitted by 
cosmological observations. Doing this one finds that 
(see also Jl6[) 



9.1 x 1CT 52 mpi = 1.1 x 10~ 23 eU : 



(2) 
(3) 
(4) 



A ~ 20.28, 
V ~ (3 x HT 27 m Pl )\ 

with mpi = G -1 / 2 the Planck mass. 

If a galaxy is an oscillaton, i.e. an oscillating soliton 
object, it must correspond to coherent scalar oscillations 
around the minimum of the scalar potential (Q). The 
scalar field $ and the metric coefficients (considering the 
spherically symmetric case) would be time dependent, 
and it has been shown that such a configuration can be 
stable, non-singular and asymptotically flat Q. For the 
scalar field collapse, the critical value for the mass of an 
oscillaton (the maximum mass for which a stable con- 
figuration exists) will depend on the mass of the boson. 
Roughly speaking, if we take m$ = 1.1 x 10~ 23 eU, and 
use the formula for the critical mass of the oscillaton cor- 
responding to a scalar field with a <I> 2 potential (i.e. just 
a mass term), we expect the critical mass to be jj],^] 



M c „ t ~0.6^~10 12 M G 

771$ 



(5) 



This is a surprising result: the critical mass of the model 
shown in is of the order of magnitude of the dark 
matter content of a standard galactic halo. 



In order to study this situation for the case of a poten- 
tial of the form (|l|), we present a numerical simulation 
of Einstein's equations in which the energy momentum 
tensor is that of a real scalar field. The scenario of galac- 
tic formation we assume is as follows: a sea of scalar 
field particles fills the Universe and forms localized pri- 
mordial fluctuations that could collapse to form stable 
objects, which we will interpret as the dark matter halos 
of galaxies. 

We evolve the spherically symmetric line element 



ds 2 = -a 2 (r, t)dt 2 + a 2 (r, t)dr 2 + r 2 dU 2 



(0) 



with a(r, t) the lapse function and a(r, t) the radial met- 
ric function. We choose the polar-areal slicing condition 
(i.e. we force the line element to have the above form at 
all times, so that the area of a sphere with r — R is al- 
ways equal to AuR 2 ). This choice of gauge will force the 
lapse function ot(r,t) to satisfy an ordinary differential 
equation in r (see below). 

The energy momentum tensor of the scalar field is 

V - - 9 -f- Q + 2V{$)\ . (7) 

We now introduce the first order variables \& = <& r and 
IT = a$. t /a. Using these new variables, the Hamiltonian 
constraint becomes 



1 — a Kgr 



[* 2 + n 2 + 2a 2 U] 



a 2r 4 

and the polar-areal slicing condition takes the form: 

a r a r a 2 — 1 2 
— — = — — | Kora V . 



(8) 



(9) 



All other components of Einstein's equations either van- 
ish, or are a consequence of the last two equations. 
The Klein-Gordon (KG) equation now reads 



a 



1 / r 2 a* 



r \ a 
a 



aa 



dV_ 



(10) 
(11) 
(12) 



Equations (J8j-y_2j) form the complete set of differential 
equations to be solved numerically. For numerical pur- 
poses, the evolution equation for II above is further trans- 
formed into the equivalent form: 



d fr 2 a^\ dV 
dr 3 V a J d$ 



(13) 



Notice that the first term on the right hand side of this 
equation includes now a first derivative with respect to 
r 3 (and not a third derivative) . The reason for doing this 
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transformation has to do with the numerical regulariza- 
tion near the origin of the 1/r 2 factor in equation (|ll| ) 
above (see Ref. p7[). 

In order to deal with non-dimensional units, we define 
x = Ir. A natural scale for the potential is given by 
l^ 1 = = \2pc = 40yr. The parameter I also 

gives the time scale r = tl. The scalar field has an initial 
Gaussian profile 



Hzq $(x, t = 0) = Ae 



2 / 2 
-X 1 8 



(14) 



with A the amplitude and s the width of the Gaus- 
sian. The physical properties describing the state of 
the system are the energy density of the scalar field 
p$ = (m^l 2 / 8ir) p s , with p s the dimcnsionlcss quantity 



1 



2a 2 

and the integrated mass 



M{x) 



21 



f Ps (X)X 2 dX. 
Jo 



(15) 



(16) 



Some details about our numerical implementation are 
in order. To integrate the KG equation numerically we 
use a method of lines with standard centered second or- 
der finite differences in space, and a third order in time 
integrator. At the outer boundary we impose a condition 
for simple outgoing radial waves. We deal with the sin- 
gularity at x = by straddling the origin and imposing 
the adequate parity conditions for each function on an 
auxiliary point at x = —Ax/2. The ordinary differential 
equations for a and a are solved using a second order 
Runge-Kutta method. The third order in time integra- 
tion has been chosen to reduce as much as possible the 
numerical dissipation in our code, which we have found 
to be crucial in order to obtain reliable results for the 
long time runs we have studied. 

The numerical simulations suggest that the critical 
mass for the case considered here, using the scalar po- 
tential (pi), is approximately jL5|Jl^] 



0.1 



"pi 



2.5 



10 13 M(T 



(17) 



The results of the numerical simulations are as follows. 
Essentially, we have found three different types of behav- 
ior for the scalar field collapse. In the first case, a generic 
feature is that scalar field distributions with an initial 
mass slightly larger than the critical mass collapse very 
violently and form a black hole. In the second type of 
behavior, fluctuations with an initial mass significantly 
smaller than the critical mass can not form stable oscil- 
latons: the scalar field is completely ejected out as the 
system evolves flgl - The third behavior corresponds to 
a case where a fraction of the initial density is spread 
out, leaving an oscillating object that appears to be sta- 
ble. This situation happens in a narrow window of initial 
conditions, between 0.05 — 1 x M crit |Q. 
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FIG. 1. Temporal evolution of the maximum value of the 
energy density. The parameters for the initial configuration 
are A — 0.01, s = 2.0 (see (|14|)), which correspond to an 
initial mass of Mi = 3.25 x 10"~M Q . See also Fig. |. 



In Fig. ^ we show the evolution in time of the maxi- 
mum value of the energy density for an initial configura- 
tion that results in the formation of an oscillaton. In 
this case, we have taken an initial configuration with 
A = 0.01, s — 2.0, which implies an initial mass of 
Mi = 3.25 x 10 12 Mq. For this run we used Ax = 0.005, 
At = Ax/10, and 10,000 grid points, which puts the 
outer boundaries at x = 50. The run was followed until 
t = 2000, some 40 light crossing times. 
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FIG. 2. Energy density of the object at t = 2000. In the 
plot we show x 2 p a (x) using a logarithmic scale. The solid line 
corresponds to a run with Ax — 0.005 and the boundaries 
at x = 50, the dashed line to a run with Ax = 0.01 and 
boundaries at x = 50, and the dotted line to a run with 
Ax = 0.01 and boundaries at x — 100. 

In Fig. H we plot the energy density x 2 p s (x) of the ob- 
ject at t = 2000, using a logarithmic scale, for three differ- 
ent runs: the run mentioned above with Ax = 0.005 and 
boundaries at x = 50, and two runs with half the resolu- 
tion (Aa; = 0.01) and with the boundaries at x = 50 and 
x = 100 respectively. Notice how the the two runs with 
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the boundaries at the same location coincide very well 
throughout the computational domain. The run with 
the boundaries twice as far agrees well with the other 
runs for x < 20, but differs significantly outside where 
x 2 p s < 10~ 6 . This indicates that at such low levels, the 
solution is dominated by boundary effects. Our bound- 
ary condition is clearly introducing numerical noise at 
such late times. 

Figure || shows the integrated mass for the initial 
and final (t — 2000) stages of the evolution. A small 
drop of ~ 0.5% in the total integrated mass can be ob- 
served, but convergence tests suggest that most (if not 
all) of this mass loss is caused by a small amount of nu- 
merical dissipation still present in our numerical method. 
This implies that the system does not radiate any signifi- 
cant amount of energy during the time of the simulation, 
which indicates that the object is very stable. 



could be a good model for the universe at cosmologi- 
cal level, here we see that the scalar field could also be a 
good candidate for the dark matter content of individual 
galaxies (as suggested in J|,|]||] ) • 
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FIG. 3. Initial (solid line) and final (dotted line) integrated 
masses. 

From the cosmological point of view, the narrow win- 
dow of initial conditions means that not all fluctuations 
will collapse into stable objects. Moreover, the collapsed 
objects will have masses of the same order of magnitude 
Affinal ~ 10 12 M Q , as it seems to be precisely the case for 
galaxies. 

Summarizing, from the results of the numerical simu- 
lations of the collapse of the real scalar field with a cosh 
potential we find many similarities with the structure 
of the halos of galaxies. The scalar field density profile 
is not singular at the center. This fact, and the values 
of the final masses obtained using the cosmological val- 
ues (|j|) and (|4|) for the parameters of the self-interaction 
potential, could correspond to objects like realistic galax- 
ies. Moreover, it is in agreement with the observational 
constraints related to the phenomenological maximum 
galactic mass pointed out by Salucci and Burkert Jl9| . 
Therefore, we expect that fluctuations of this scalar field, 
due to Jeans instabilities, will in general collapse to form 
objects of the order of the mass of the halo of a typical 
galaxy. 

We have shown before ||,[?J that the SFDM model 
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